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Abstract 


A propeller lifting-surface design and analysis program is improved upon by implementing 
enhancements in the source distribution calculation to represent the blade thickness. It is 
recognized that the present method of setting the source line distribution representing blade 
thickness (currently based on linearized slender-body theory for an isolated foil section) may 
introduce significant errors. This is the case for propulsors with a combination of a large 
thickness/chord ratio (blockage effect) and numerous blades (cascade effect). 

A source panel (area) method was developed to more accurately model these effects. This 
method uses the lattice structure of the current PBD-14 code from which to compute the source- 
induced velocity factors between the blades, hub, and duct, if present. Using the method 
of images allows the hub and duct to be modeled as panel images from the blade panels. 
The source-induced effects of the whole propulsor are accounted for by using a panel method 
to obtain a source distribution along the mean camber surface of the blade. Invoking the 
kinematic boundary condition on the true blade suction and pressure surfaces solves this system 
of linear equations, which represent the blade thickness distribution. This robust formulation 
assigns source strengths more accurately over a much larger range of thickness/chord ratios 
and increasing numbers of blades, as evidenced by a more accurate velocity streamline trace 
representation of the actual pressure and suction side surfaces of the blade. Experimental 
validation is demonstrated for open and ducted flow stators. 
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Nomenclature 


Mathematical Notation 


A area of quadrilateral panel (flat z — plane) 

BB blade-to-blade influence functions 

DB duct-to-blade influence functions 

dR radial difference between source element endpoints 

ds differential component of local coordinate position 

dS differential component on the blade surface 

dx axial distance between control points in a chordwise strip 
dt ot differential thickness between control points in a chordwise strip 
ductrat Z aes duct element imaged source strength adjustment factor 
D propeller diameter 

hy, he ox metric coefficient 

HB hub-to-blade influence functions 

AIF horseshoe influence function matrix 

hubrat pildelem. hub element imaged source strength adjustment factor 
L pressure side of the blade (lower surface) 

Lotdelem length of blade source element 

ceductelem length of blade duct element 

Lhubelem length of blade hub element 

Ltotat total length of a source element 

m mean camber surface of the blade 

nH normal vector on a surface 


Pp surface coordinate 
P arbitrary point in space 
R radial coordinate 


X(S1, 82) — X(s{,85) differential position vector 


81,82 local coordinates 
S blades surface 
SIF source influence function matrix 
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oe 


wake surface created by horseshoe 
tangential vector on a surface 

velocity associated with the x-coordinate 
inflow velocity to the propeller - 
component of inflow velocity 

velocity associated with the y-coordinate 
effective velocity 

total induced velocity 

induced circumferential mean velocity 
fluctuating component of induced velocity 
total velocity 

component of velocity due to vortices 
component of velocity due to sources 
velocity associated with the z-coordinate 
surface coordinate vector for x, y, z 

axial coordinate 

vertical coordinate 

z-coordinate 

upper/lower blade surface (suction/pressure side) 
angle of attack of propeller blade 
potential 

vortex element strength (per unit length) 
blade circulation 

blade rotation rate (224) 

source panel distribution (per unit area) 
source strength distribution (2D sense) 
source line distribution (per unit length) 
sum of positive source strengths 

sum of negative source strengths (sinks) 
sum of panel areas due to source strengths 


sum of panel areas due to sink strengths 
11 


Abbreviations 

CMV circumferential mean velocity 

horseshoe a vortex structure composed of: a spanwise vortex, 
its shed chordwise vortices on the blade and in the 
transition wake, and the ultimate wake helices. 


RANS Reynolds-Averaged Navier Stokes 
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Chapter 1 


Introduction 


1.1 Overview 


The goal of this research was to create a more accurate model which describes the true lifting 
surface of a propeller blade system. This system includes the interactions between the propeller 
blades and the hub, as well as the duct, if present. Since many current propeller system design 
trends include the use of a multiple-staged ducted propulsors, accurately modelling the effects of 
induced velocities between the different components of the propulsor are paramount to perform 
a good analysis. There are several numerical tools and techniques which accomplish this task 


sufficiently [1]: 
1. Lifting-Line Methods 
2. Vortex-Lattice Lifting-Surface Methods 
3. Boundary Integral Equation Methods 


4, Axisymmetric and Three Dimensional Reynolds Averaged Navier Stokes (RANS) Flow 


Solvers. 


The vortex-lattice lifting-surface technique provides an accurate and robust method to examine 
the design [2] and analysis [3] of complex marine propulsors efficiently. Kerwin et al. [2] 
[3] developed such a lifting-surface technique which has been improved upon [4] [5] since its 


introduction. 
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The vortex-lattice lifting-surface technique is a potential flow method which does not fully 
describe what occurs within the propulsor system, like full viscous flow models. However, 
Kerwin et al. [4] developed a coupling procedure which uses the potential flow lifting-surface 
technique in conjunction with a RANS viscous flow solver. This coupled technique allows for 
multiple blade row analysis while providing an accurate model of the physical behavior of the 
system. 

The linearized slender-body theory representation does not accurately model the flow over 
the propeller blades because it under predicts the effects of the thickness distribution of the 
propeller blades. It also does not account for the blade-to-blade, blade-to-hub, and blade-to- 
duct interactions. Therefore, as more blades and stages are added to the propeller system, 
the linearized slender-body theory method is unable to accurately predict the cascade and 
blockage effects created. This deficiency in the current vortex-lattice method is the motivation 
to implement an enhanced thickness model in an effort to better approximate the flow patterns 


through the system. 


1.2 Background 


1.2.1 Vortex-Lattice Lifting-Surface Analysis Tools 


A comprehensive explanation of the vortex-lattice lifting-surface theory and methods for the 
design and analysis of marine propulsors is documented extensively (2], [3], [6], and [7]. Propeller 
Blade Design (hereafter designated PBD-14) developed by Kerwin! et al. was used for this 
research. In PBD-14, a lattice geometry of discrete vortex segments is placed along the blade 
mean camber surface. Control points are placed within each grid lattice. Blade thickness is 
represented using a linearize slender-body theory to obtain a source strength distribution at 
the lattice segments in a strip-wise sense. Each lattice segment has known source strength 
and the vorticity strength is solved for, such that the kinematic boundary condition is satisfied 
along the mean camber surface, namely the total velocity must be tangent to the surface at 
each control point. A trailing wake system is grown from the trailing edge of the blades. See 


Figure 1-1 for the vortex-lattice discretization of the propeller geometry. 


‘Justin E. Kerwin, Profesor of Naval Architecture, Ocean Engineering Department, MIT. 
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ieee 
,|Blade Vortex Lattice |, 
with Control Points 


Hub Vortex Transition Wake 
Lattice 


Figure 1-1: Vortex lattice discretization of a propeller geometry (4119) where the symbols 
represent control points. , 


1.2.2 Propeller Blade Thickness Modeling 


It is recognized that the present method of setting the strength of the source distribution 
representing the blade thickness (based on linearized slender-body theory for an isolated foil) 
may introduce significant errors in the case of propulsors with a combination of a large thick- 


ness/chord ratio and numerous blades. The questions to be answered are: 


1. Is the thickness-induced velocity correct in PBD-14? There are three parts to this prob- 


lem: 


e isolated two-dimensional foil, 
e two-dimensional cascade effects (emulations of multi-blade flow) and 


e three-dimensional effects. 
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2. How does one accurately model these effects in a marine propulsor system? Two methods 


will be discussed below: 


e linearized slender-body theory and 


e panel method for a source distribution along the mean camber surface. 


Due to efforts to more accurately assign source strengths to the vortex-lattice structure 


(i.e. correctly model cascade and blockage effects in propulsors), enhancements were sought in 


the current methodology of the source distribution to represent the blade thickness distribu- 


tion. The current version of PBD-14 uses linearized slender-body theory to predict the source 


strength necessary to satisfy a pseudo kinematic boundary condition at the designated control 


points. Regardless of the blade geometry, the kinematic boundary condition is invoked for flow 


over a thickened foil, assuming zero camber and loading. It performs a systematic assignment 


of source strengths, one chordwise strip at a time using equations 1.1 and 1.2: 


2 2 St 
oL= / odz = Viotal aint (1.1) 
Li xy ah 


which simplifies to 


where 
o 
OL 
Viotal 
ét 
dx 


Ill 


ill 


ill 


OL = Vota dt (1.2) 


the source strength distribution 

the source strength distribution per unit length 

the total inflow velocity in the presence of an operating propeller 
the increase in blade thickness between lattice control points 
along a given chordwise strip 

the axial distance between control points being compared to 


determine the source strength 
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0 0.025 0.05 beg 01 0.125 0.15 


Figure 1-2: Comparison of desired thickness form and source-induced thickness form of a Bi- 
convex blade section obtained by streamline tracing for current PBD-14 code. 


Therefore, the strength of each source line is determined using the relationship between the 
total inflow velocity and the differential thickness of the blade at some position (x) downstream 
from the leading edge of the blade section. Figure 1-2 shows an extreme test case for an infinite 
pitch stator in uniform flow with zero circulation (I) and a thickness/chord ratio of 27 percent. 
The desired thickness form (solid line) can be compared with the actual thickness form (vector 
symbols) obtained by streamline tracing. The source-induced thickness form of the section is 


about 20 percent too thin, thus indicating the source strength is too small. 


Source Line Method 


In this method, the source distribution is determined for the vortex-lattice structure by modeling 
the complete blade system (all blades, the hub, and the duct, if present). This is a vortex-lattice 
method which distributes the source line elements on the mean camber surface of a propeller 
blade. Then, the source distribution for the entire blade system is determined by satisfying 
the kinematic boundary condition at control points which are on the actual blade surface 
(pressure and suction sides). The result is a source strength per unit length at the vortex/source 


elements. This yields extremely accurate results for the foil shape shown in Figure 1-2, but gets 
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progressively inaccurate once the thickness/chord ratio falls below ten percent. 


Source Panel Method 


A alternative and more robust method, determines the source distribution for the complete 
blade grid structure by fully modeling the entire blade system (all blades, the hub, and the 
duct, if present) with source panels. This is a panel method which distributes the source panels 
on the mean camber surface of a propeller blade. Then, the source distribution for the entire 
blade system is determined by satisfying the kinematic boundary condition at control points 
which are on the actual blade surface (pressure and suction sides). The result is a source 
strength per unit area which can be integrated over the blade grid structure to obtain a similar 
lattice source strength per unit length as above. This yields extremely accurate results for the 


foil shape shown in Figure 1-2 for all ranges of the thickness/chord ratio. 


1.2.3. Hub and Duct Modeling 


The hub is represented by a similar vortex-lattice geometry which is appropriately placed using 
the method of images. The duct, if present, is developed in the same manner. The constant- 
strength source line elements in the hub and duct are ratioed from the blade source element 
strengths by comparing the lengths of the blade elements with their associated imaged elements 


in the hub and duct. Figure 1-3 shows a hub and duct representation in PBD-14. 


1.3. Objective 


This thesis provides a description of the advanced thickness modeling techniques implemented 
into the lifting-surface propeller blade design and analysis code PBD14.8 [8]. This enhancement 
extends the capability of current lifting-surface theory by properly accounting for thickness 
effects (cascading and blockage) in this model structure more completely. An overview of the 
process is presented first, followed by more detailed descriptions of the thickness modeling 


process. Finally, validation examples are presented. 
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Blade Lattice 


HH 


Hub Image 
Lattice 


Figure 1-3: Blade, hub and duct representation in PBD-14 for Ka-455 propulsor. 
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Chapter 2 


Lifting-Surface Design and Analysis 
Theory 


2.1 Overview 


In the past, lifting-surface theory has provided accurate results by using slender-body theory 
thereby approximating higher-order thickness and higher-order viscous effects for conventional 
propulsors as equivalent. It has been shown that in conventional thickness and loading distri- 
butions, the decrease in lift due to viscous effects is offset by the increase in lift due to blade 
thickness [2]. Due to the advances in blade section design, the higher-order viscous and higher- 
order thickness effects do not necessarily cancel. Hence, these factors must be accounted for in 
the design and analysis of marine propulsors. 

Lifting-surface propeller design and analysis places a vortex-lattice structure on the mean 
camber surface of the blade. The presence of thickness is accounted for by the distribution of 
sources and sink lines collocated with the vortex elements. There is a local force as a result 
of the velocities of a source, influencing its neighboring vortex elements. An inter-blade effect 
occurs due to the blockage of the modeled thickness of the blade. There is also a Lagally force 
that arises from the force on the sources themselves. Panel method codes model the thickness 
of the blade, so these effects are addressed. However, current lifting-surface codes do not have 
this ability [3] [9]. Also, blade thickness results in a blockage to the axisymmetric flow which 


affects the flow distribution past the propulsor and duct, if present. This factor must also be 
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taken into consideration if propulsor performance is to be predicted more accurately [10]. 


2.2 Vortex-Lattice Method 


The vortex-lattice method of solving the lifting-surface propeller blade problem can be catego- 


rized into four separate steps: 


1. A continuous vortex/source sheet, located on the blade’s mean camber surface, is linearly 


represented by a discretized geometry for the propeller blade. 
2. Determine the influence coefficients for the vortex/source distributions. 


3. Apply the kinematic boundary condition which stipulates that the flow on the blade 


surface is only in the tangential direction. 


4. Solve the system of linear equations developed from the steps outlined above. 


2.2.1 Geometry Discretization 


Numerically, circulation and thickness distributions are discretized using a lattice structure 
of constant-strength, straight-line vortex/source elements that are aligned in the chordwise 
and spanwise directions. These vortex/source elements are collocated along the blade’s mean 
camber surface. PBD-14 uses a vortex-lattice technique which incorporates a uniform B-spline 
representation of the blade’s mean camber surface [2]. This representation of the blade allows 
for efficient interrogation of the B-spline control polygon to discretize the blade’s surface [6]. 
Furthermore, various lattice structure spacing options as available (i.e. cosine and/uniform 
spacing) as described by Greeley and Kerwin [3]. Kelvin’s Theorem states that vortex lines 
cannot end anywhere in the flow field, therefore trailing vortices convect downstream. The 
bound element with its associated convected trailers form a horseshoe vortex which extends 
beyond the trailing edge of the propeller blade [3]. The trailing vortices must be aligned with 


the flow in order to be force-free. 


21 


2.2.2 Influence Coefficients 


Lifting-surface methodology solves the propulsor problem by first developing an influence ma- 
trix, [HIF; ;]. This matrix uses the induced velocity produced by each vortex/source-lattice 
segment as the influence (assumed unit strength of one) upon each of the blade control points. 
Therefore, the influence coefficients depend solely on the geometry. When this influence matrix 
is multiplied by the actual vortex segment strength, the induced velocity at every control point is 
known. Section 2.2.4 more fully explains the influence coefficients purpose in the lifting-surface 


methodology. 


2.2.3. Kinematic Boundary Condition 


In order to correctly describe the flow through the propeller blade system, one must implement 
the appropriate boundary condition on the rigid body. The kinematic boundary condition 
stipulates that no flow will penetrate the surface of the blade (ie. the total velocity flow is 
tangential at the surface). Another way of stating the kinematic boundary condition is that the 
total velocity normal to the grid surface must be zero. This condition is specified at a number 
of control points on the mean camber surface which is equal to the number of vortex/source 


elements representing the blade lattice structure and is shown in equation 2.1 below: 


Viotal -n=0. (2.1) 


Section 2.2.4 more fully explains how the kinematic boundary condition is utilized to solve the 


lifting-surface problem. 


2.2.4 Solution Procedure 


The lifting surface can be defined by a curvilinear system of coordinates (51,52) on the blade 


surface as shown in Figure 2-1. The surface coordinates are defined as 


x(1, $2) = [x(s1, 82), y(S1, $2), 2(81, $2)] (2.2) 


and the metric coefficients are given as 
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Ox 
epee 2.3 
acre (2.3) 
and 
Ox 
eR eS 2.4 
eal (2.4) 
‘ or Sy 
,|Blade Vortex Lattice |; 
with Control Points 
Hub Vortex Transition Wake 


Lattice 


Figure 2-1: Vortex lattice discretization of a propeller geometry where the symbols represent 
control points and s; and sg are orthogonal vectors on the blade surface. 


The kinematic boundary condition is given as 


81,85) xX R|-n ee 
| = —[Uin(s1, $2) + OR + Vina(S1, $2)] - fi (2.5) 


where 
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R = x(s1,82) —x(s{,8%) 


=> 
ll 


vector normal to S at point (sj, 89). 


Equation 2.5 specifies that there is no flow through the blade surface S, considering a 
continuous sheet of vorticity (y) on the blade surface and in the wake (S,,) in the presence of 
an incoming flow (Ujn) and rotational velocity (QR). Since the thickness distribution for the 
blade is known, the source-induced velocities (Vind,,,..) may be computed and included on the 
right hand side of the equation [10]. 

A lifting-surface code is used for either a design of a new propulsor or an analysis of an 
existing propulsor. In each case, equation 2.5 applies, but the known versus unknowns vary. 
For propulsor design, a desired radial and chordwise loading distribution is prescribed. Then, 
the blade shape is manipulated until equation 2.5 is satisfied. For propulsor analysis, the 
blade shape is prescribed. Then, the resulting circulation is solved. This process summation is 


simplified to matrix equation 2.6: 


(HIF; ;| . [P5] =, —[V; ‘nl (2.6) 
where 

HIF; = matrix of horseshoe influence coefficients. The subscripts denote the influence of 
the j*" horseshoe vortex on the i** control point. Each element represents the 
velocity component normal to S, induced at a control point due to a horseshoe 
vortex of unit strength. 

[T; = blade circulation scalar composed of individual vortex strengths of the jh 

vortex. 

V;-fi = velocity component normal to the blade mean camber surface at the i” control 


point. 
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Solving this matrix system of linear equations yields equation 2.7: 


[Ps] = [Vj AIF, J (2.7) 


Blade thickness effects are added by placing discrete source lines coincident with the chord- 
wise blade vortex elements. The resulting propulsor forces due to the lifting surface and thick- 
ness are calculated from the Kutta-Joukowski and Lagally theorems, respectively. A leading 
edge suction force and Lighthill pressure distribution correction are applied to those forces, and 
the propulsor’s sectional viscous drag is calculated based on either stripwise two-dimensional 


empirical drag coefficients or a stripwise two-dimensional integral boundary layer calculation. 


2.3. Propulsor Inflow 


The accuracy of modeling the propulsor is dependent upon correctly representing the propeller 


inflow. The three concepts for propeller inflow are: 


Definition 1 Nominal inflow: Velocity field at the plane of the propeller when the propeller is 


not operating. 


Definition 2 Total inflow: Velocity field at the plane of the propeller when the propeller is 


operating. 


Definition 3 Effective inflow: Total inflow without the effects of the propeller induced veloci- 


ties. 


In propeller design and analysis, the lifting-surface methodology deals with the flow resulting 
from the operation of the propeller in a specified axisymmetric effective inflow. Therefore, the 
total velocity relation at the blade surface in the axial, radial and tangential directions can be 


decomposed into the following three components 


Viotal ae Veff ae Virgen + V Mac (2.8) 


where 
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Viotal 


Veff 
Vind 


CMV 


Vindptuc 


Ill 


the total velocity in the presence of an operating propeller at a given point on 
the blade. 

the effective wake behind the ship. 

the circumferential mean velocity induced due to propeller blades and wakes. 


the fluctuating component of induced velocity due to propeller blades. 


Therefore, the total velocity is a combination of the effective inflow and the net propeller 


blade induced velocities. In the propeller blade problem, the total velocity is used to solve the 


kinematic boundary condition as shown in equation 2.1. 
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Chapter 3 


Numerical Implementation of the 


Enhanced Blade Thickness Model 


3.1 Overview 


In the current version of PBD-14, propeller lifting-surface design and analysis program code, 
the thickness distribution is developed in the following manner using linearized slender-body 


theory: 


1. For a chordwise strip, determine the thickness distribution from the B-splined input file 
parameters and the change in thickness (6¢) between the vortex-lattice structure control 


points (dz). 


2. Use the flow field’s circumferential mean velocity (CMV) and rotational velocity compo- 
nent ((Q2R) to determine the blade’s local velocity at the control points (along the mean 


camber surface). 


3. Solve for the source/sink line elements constant-strength (o7;,;) necessary to create the 
correct 6¢t between the control points on that spanwise strip. Equation 3.1 shows that the 
source strength per unit length is a function of inflow velocity, the change in the thickness 


distribution and the differential length between control points: 


OLijg = Viotal dt. 
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4. Adjust each source line element strength for “slant.” This effect increases the strength of 
the source line element when it is not perpendicularly aligned with the inflow by a factor 
of Etta Here, Lictat is the total length of the element and dR is the radial difference 


between element endpoints. 


5. Correct the sink strengths to ensure each chordwise strip of source/sink line elements 
sums to zero. This effect increases each sink element strength by a factor of ane thus 


ensuring that the 3° Ogre + > Osnz = 0. 


This methodology was flawed for cases where the blade design had numerous blades and/or 
a high thickness/chord ratio (> 10%). This scheme did not account for blade-on-blade, hub- 
on-blade, or duct-on-blade effects in the source line strength determination. Therefore, as the 
source lines of each modeled propeller section and its associated components get closer together, 
these effects increase in importance (i.e. far field effects begin to move into the near field range 
to increase the influence of these elements on one another). 

Because of the shortcomings of this methodology, more accurate methods of accounting for 
blade, hub and duct interactions to more fully represent the propulsor were investigated. Two 


promising techniques were developed, implemented and tested in PBD-1/: 
1. Source Line Method 
2. Source Panel Method. 


The source panel method is an improvement upon the work completed in the source line 
method. This work was undertaken to develop a more robust representation of all the effects 
that determine the blade thickness distribution. For example, the source line method accounted 
for the cascade effect (increasing the number of blades creates larger blade-to-blade influences). 
It also was able to properly represent the blade thickness distribution for large thickness/chord 
ratios. Both of these efforts were improvements over the shortcomings of the previous PBD-14 
methodology of assigning source strengths. However, the source line method over predicted the 
blade thickness distribution when the thickness/chord ratio decreased to less than ten percent. 
This effect got worse as the thickness/chord ratio decreased. The source panel method corrects 
the shortcomings of the source line method, so that it accurately represents the blade thickness 


distribution over the broadest range of thickness/chord ratios and numbers of blades. 
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3.2 Source Line Method 


This thesis seeks to find improvements in PBD-14 to more accurately represent the blade 
thickness distribution over the blade lattice structure in the presence of numerous blades or 
large thickness/chord ratios led to the development of the source line methodology. An outline 


of this method follows: 


1. Discretize the blade geometry on the mean camber surface so that the vortex and source 


elements are collocated. 


2. Develop each chordwise blade thickness distribution from the PBD-14 input file 4 values 
shown in Appendix B. 


3. Fully develop the blade’s actual outer surface for each spanwise position at each vortex 
element endpoint that lies on the mean camber surface. This is accomplished by creating 
the pressure and suction side surfaces from the actual thickness distribution of the blade 
section at the mean the camber surface vortex/source elements. The mean camber surface 
lattice endpoints were moved in the normal direction by the appropriate 6t for the specified 
point. Figures 3-1 and 3-2 show the discretized geometry of a propeller blade with the 


pressure and suction side surfaces included. 


4. Develop a complete set of control points and a normal to the control points on the pressure 
and suction side of the actual blade surface. This is done in the same fashion as the mean 


camber surface calculations in previous versions of PBD-14. 


5. Develop the hub and duct imaged lattice structure. The hub and duct images lie along 
the mean camber surface of the inner and outermost spanwise set of endpoints, thus 
having the same pitch of the blade at the hub and duct intersection points. Also, the 
hub and duct imaged source strengths are adjusted by a factor of Ftldelem. and —Ptldelem 


Lnubelem Lauctelem’ 


respectively. 


6. Determine the total inflow velocity at the control points on the blade’s mean camber 


surface. 
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7. Develop a system of linear equations in matrix form, similar to equation 2.6, to solve 
for the source strengths associated with each of the source-induced function coefficients 
[SIF;,;|. This is done while satisfying the kinematic boundary condition (Victa - i = 0) 
at the control points on the actual pressure and suction side surfaces of the blade. The set 
of matrix equations is increased to add a further condition that the 5° Osre + > Osnk = 0 


for each chordwise strip on the blade. 


This set of matrix equations is shown in a simplified notation in equation 3.2: 


[SIF; ] ‘[o1,) =—[V;- a (3.2) 

where 

SIF;; = matrix of source influence coefficients. The subscripts denote the influence of 
the j* source element on the i*” control point. Each element represents the 
velocity component normal to S, induced at a control point due to a source 
element of unit strength. 

o,, = constant strength source composed of individual source strengths of the hg 

source elemnts. 

V;-fi = velocity component normal to the blade mean camber surface at the 7” 


control point. 


Solving this matrix system of linear equations yields equation 3.3 


(o3] = —[V;-a)[SEF, (3.3) 


In order to determine what truly affects the blade thickness distribution problem, a decom- 
position of source and vortex velocities at the pressure, suction and mean camber surfaces was 
necessary. Figure 3-3 shows a simple chordwise blade representation for the velocity effects on 


the blade thickness distribution. The proof follows. 
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Figure 3-1: XY-coordinates for blade thickness discretization on the pressure and suction sides. 


Definitions for the nomenclature in Figure 3-3 is: 


— = pressure surface (lower) of the blade. 


+ = suction surface (upper) of the blade. 
m = mean camber surface (middle) of the blade. 
n = normal for pressure and suction side surfaces. 
t = tangent for pressure and suction side surfaces. 
(s1,82) = local coordinates in normal (1) and tangential (2) 
directions on the mean camber surface. 
(U;,,Us,) = component of inflow velocity in the normal (1) 


and tangential (2) directions. 


ill 


(Vos, ? Vo5q) 


Vans VV) 


component of velocity due to sources. 


component of velocity due to vortices. 
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Figure 3-2: X Z-coordinates for blade thickness discretization on the pressure and suction sides. 


In order to meet the kinematic boundary condition on both the pressure and suction side 


surfaces of the blade, equations 3.4 and 3.5 must hold true. 


[War + Vor + V 54) H+ (Usp + Voy + Ve) =0 (3.4) 
[(Us, + Von, + Ve) + Use + Von +V%e,)t]_ = 0 (3.5) 
where 
fi, = —fi_ and ty = t_ exactly, 


and for small thickness distributions, 
Va & V¥51_. 
From substitution and adding equations 3.4 and 3.5, the following result is obtained 


(Vos, 5 Vos.) ae 


Von, +Vem) , WV, +n) $=0 (6) 
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Figure 3-3: Chordwise strip of a blade section to represent the components which affect blade 
thickness distribution. 


Assuming a small thickness distribution, Figure 3-4 shows a simplified mean camber line rep- 
resentation of V7,,. 
Since Vy, =V%s.,, + and V¥s0e = V Vso m — 4) then averaging the effects of the vortex 


induced velocities from the pressure and suction surfaces of the blade at the mean camber 


V¥5q, tV 7 
surface produces Vy,, = a = V7,, ,: Hence, equation 3.6 can be simplified further 


(Vos, —Vos_). 


2 ny + [bs + Worn, + Vor_) 


; +V7,q,{ t4 =0. (3.7) 

Furthermore, for a symmetric, two dimensional foil, V7. = 0. Also, if the foil has a moderate 
camber, like most propeller blade sections, V5. * 0 and can be considered as a part of the 
CMV inflow. This assumes that V¥,, is a small percentage of total CMV inflow which is made 


up of many components. Therefore, the Vs, component of velocity can be moved to the right 
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Figure 3-4: Mean camber line representation of vortex-induced velocities averaging to Vy,,. 


hand side of equation 2.6 and iterated upon to get a final solution. 

Finally, Vos,, and Vos,_ on the pressure and suction side surfaces of the blade can be 
expressed as a dot product of the source-induced influence coefficients [SIF ;] and the unknown 
source line strength [o,,|. Thus, equation 3.7 simplifies to the matrix equation 3.2. 

To check this reasoning, a symmetrical, two-dimensional foil is placed at an arbitrary angle 
of attack (a) as shown in Figure 3-5. Since Us, , = Uinsina, Us, = Uinsina and Us, = 


Us,_ = Uin cosa, then substituting these parameters into equation 3.7, yields the expected 


result that the blade thickness distribution is only affected by Ujn cosa or 


Vos, fit + Uin cosa + Voz.) ty =0. (3.8) 


Once source-induced coefficients are developed for the pressure and suction side surfaces 


of the blade, they are averaged to get a representative [SIF; ;] coefficient matrix for the mean 
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a I, o 
S> i 


Figure 3-5: Two-dimensional blade section used to check what attributes affect blade thickness 
distribution. 


camber surface. The mathematical formulation is solved using a least squares methodology 
which already existed in the PBD-14 code. In this solution formulation, the hub-to-blade and 
duct-to-blade imaged influence coefficients are added directly to the blade-to-blade influence 
coefficients. This is explained more fully in Section 3.4.2. This is a more robust and accurate 
method than the current PBD-14 code since the kinematic boundary condition is satisfied on 


the actual blade surface at the new set of control points developed on this surface. 


3.3 Source Panel Method 


A second method to accurately model the source distributions over the blade’s mean camber 
surface was developed after the shortcomings of the source line method were found. While 


the source line method accurately represents propulsors with numerous blades and large thick- 
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ness/chord ratios, it fails to correctly represent thinner blades where the thickness/chord ratio 
approaches zero. In fact, it begins to over predict the source line strengths and hence the blade 
thickness distribution when the thickness/chord ratio becomes less than about ten percent. 
Because source line elements are a set of discrete concentrations of constant-strength sources, 
they tend to average the effects between control point positions, where the kinematic boundary 
condition is invoked for equation 3.9. This becomes more difficult for the system of linear 
equations as the blade thickness distribution approaches zero (i.e. matrix equation 3.2 becomes 
ill-conditioned and obtaining a solution by the least squares method is impossible). Therefore, 
the next logical step to obtain a less discrete concentration of the constant strength source 
distributions, but still accurately represent the blade thickness distribution is to move from line 
sources (a7) to area sources (a4) or panels. This source panel methodology seeks to obtain 
results with the same accuracy as the source line method at large thickness/chord ratios and 
numerous blades. However, it also seeks to improve upon the linearized slender-body theory 
approximation results of the current PBD-14 code as thickness/chord ratio approaches zero. 


This source panel method is outlined below: 


1. Discretize the entire B-spline blade geometry on the mean camber surface, to include both 
the leading and trailing edges. This means that one extra control point must be created 


at the leading edge. 


2. Develop each chordwise blade thickness distribution from the PBD-14 input file 5 values 
shown in Appendix B. 


3. Fully develop the blade’s actual outer surface for each spanwise position at each vortex 
element endpoint that lies on the mean camber surface. This is accomplished by creating 
the pressure and suction side surfaces from the actual thickness distribution of the blade 
section at the mean the camber surface vortex/source elements. The mean camber surface 
lattice endpoints were moved in the normal direction by the appropriate 6¢ for the specified 
point. Figures 3-1 and 3-2 show the discretized geometry of a propeller blade with the 


pressure and suction side surfaces included. 


4. Develop a complete set of control points and a normal to the control points on the pressure 


and suction side of the actual blade surface. This is done in the same fashion as the mean 
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camber surface calculations in previous versions of PBD-14. However, the leading and 
trailing edge control points are moved off the actual edge of the mean camber surface to 


a position between the lattice structure elements. 


5. Develop the hub and duct imaged panel structure. The hub and duct images lie along the 
mean camber surface of the inner and outermost spanwise set of endpoints, thus having 
the same pitch of the blade at the hub and duct intersection points. However, the hub and 
duct imaged source strengths no longer require an adjustment by a factor of Ltidelem and 


Lnubelem 


fiddle respectively. This is due to the robust nature of the panel methodology imposed 


with this technique. 


6. Determine the total inflow velocity at the control points on the blade’s mean camber 


surface. 


7. Develop a system of linear equations in matrix form, similar to equation 3.2, to solve 
for the source strengths associated with each of the source-induced function coefficients 
[SIF;,;]. This is done while satisfying the kinematic boundary condition (Viotai - fi = 0) 
at the control points on the actual pressure and suction side surfaces of the blade. The set 
of matrix equations is increased to add a further condition that the 37 Asre + 32 Asng = 0 


for each chordwise strip. 


This set of matrix equations is shown in a simplified notation in equation 3.9: 


[SIF; ;] -[o4;] = —[Vj - Al. (3.9) 


Similar to the reasoning discussed in Section 3.2, the same conclusions can be drawn about 
the velocity contributions to the blade thickness distribution. The following section will cover 
the modifications required to the PBD-1{ code to implement this methodology. 

New influence function subroutines have to be coded with the equations found in Appendix 


A. In the three-dimensional case, the discretization has two parts: 


® discretize the geometry 


e determine the singularity element distribution. 
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If these elements are represented by a B-spline net, both geometry and singularity strength, 
then the first order approximation to the surface can be defined as a quadrilateral panel! 
with a constant-strength source singularity. Since a vortex-lattice mean camber surface is 
represented as a lattice structure of rectilinear panels,? a conversion from a rectilinear surface 
to a quadrilateral surface is performed [11]. From Appendix A, the potential at an arbitrary 


point P(«,y,z) due to this quadrilateral element is 


OA dS 
Oz, y, Z = 8 | aes 3.10 
MES aI (x — 29)? + (y ~— yo)* + (z — 20)* oe 
and the velocity components can be obtained by differentiating the velocity potential: 
OD OD OF 
(u,v,w) = (57, Oy’ Be? (3.11) 


This differentiation from equation 3.11 results in source-induced velocity coefficients used to 


obtain blade-to-blade, hub-to-blade, and duct-to-blade influence functions. 


3.4 Integration of the Source Panel Model in PBD-14 


3.4.1 Blade Lattice Modification 


The current methodology uses a vortex-lattice structure which does not include the actual 
leading edge section (B-spline leading edge endpoints to the first vortex/source element). The 
new methodology is to use the entire cubic B-spline blade structure, leading edge to trailing 
edge, to describe the source distribution lines. This allowed for a symmetric distribution of 
source panels for cosine spacing in the chordwise direction along the mean camber surface. 
However, the number of source panels increases by one more than vortex elements. Control 


points were placed as before with two exceptions: 


1. a leading edge control point was added and placed in a similar fashion as the other control 


points (about the middle of the lattice rectilinear panel); and 


‘A quadrilateral panel is a flat surface with four straight sides. 
“A rectilinear panel has straight but not necessarily flat sides that can be twisted. 


2. the trailing edge control point was moved and placed in a similar fashion of the other 


control points, vice right on the trailing edge. 


These two exceptions were necessary so that the control points did not lie on the edge of a 
source panel and create a singularity for the system of linear equations, thus forming an ill- 
conditioned matrix. Similarly, the hub and duct images were developed from source panels vice 
source line elements. 

To further enhance the model and truly satisfy the kinematic boundary condition, the lattice 
points for the entire mean camber surface were moved to both the pressure and suction side 
surfaces of the blade. From these new surfaces, control points were added in the same fashion 
as the current PBD-14 code did on the mean camber surface. Figures 3-1 and 3-2 show the 
fully developed blade lattice structure. 

Lastly, the rectilinear panels which form the entire mean camber surface set of source panels 
had to be transformed into quadrilateral panels which approximate the same constant-strength 
source distribution. This transformation was completed so that the code in reference [11], 
which solves for 04, could be utilized. This transformation was accomplished using a set of 
subroutines for PSF10.3 code. These subroutines input the rectilinear panel ends consisting 
of four endpoint coordinates. Then, the subroutines return the quadrilateral panel endpoints 
which lie on a constant z — plane and a set of coordinates for the actual center of the original 
rectilinear panel. Therefore, the orientation between the local and global coordinate system is 


maintained. 


3.4.2 Adjustments to the Solution Procedure 


The solution procedure for constant-strength source area distribution is completely different 
from the current PBD-14 code. It now accounts for blade-to-blade, hub-to-blade and duct- 
to-blade influence interactions which affects induced velocities produced within the propulsor. 
These new influence coefficients are the averaged values of the pressure and suction side [SIF;,;] 
matrix which represents the mean camber surface [SIF,,;] matrix. The solution is solved as a 
system of linear equations via a least squares method vice a strict, two-dimensional, linearized, 
slender-body theory. By incorporating the effects of the other blades, the hub and the duct, 


if present, the model more accurately represents cascade and blockage effects which affect the 
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blade thickness distribution. Therefore, the propulsor blade representation is more accurate in 
a larger variety of propulsor types. 

Since panels are now utilized to obtain the constant-strength source distributions, there 
is a discrepancy in how PBD-14 must account for this source distribution. One method is 
to integrate the source panel strengths in the chordwise direction. This integration technique 
would lump the effects of the panels into source line elements which are collocated with the 
vortex elements. While this method is still an improvement over the source line method and 
linearized slender-body theory, it still diminishes the robust nature of the method. 

Therefore, a separate velocity matrix was developed from the new known source strength 
distribution for the blade, hub and duct interactions at the old control point positions on the 


mean camber surface. 


[Vind; fi] = (SIF, ,] -[o4,]- (3.12) 


Equation 3.12 shows the velocity matrix that will be imported into the [HIF;,;] system of 
linear equations, as a known velocity vector of source-induced effects on the right-hand side of 
the equation. This velocity vector will be calculated using the source panel influences on the 
specified control points vice the source line element influences. Then, this modified version of 


the [HIF;,;] system of equations will be solved to obtain [Ij]. 


3.4.3 Adjustments to Hub and Duct Modeling 


The method of images is used in PBD-14 to represent the hub and duct. This provides a means 
to implement boundary conditions necessary in potential flow theory. For singularities on rigid 
boundaries, potential flow introduces another singularity into the flow field which mirrors the 
original singularity. The vortex element images model the hub and duct accurately and are 
satisfactory from the standpoint of meeting the kinematic boundary condition. However, the 
source element images are less accurate. The resulting imaged-lattice structure is exhibited in 
Figure 1-3. The imaged source panels for the hub and duct are created in the same manner as 


the blade source panels. 
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Role of the Hub and Duct 


Adding the hub and duct source influence functions to the [SIF;,;] matrix system expands 
the complexity of appropriately accounting the influences on each control point, but does not 
complicate the actual system of linear equations that must be solved. Equation 3.13 presents 


the expanded influence function coefficient for a generalized control point position: 


Source Panels 
Control | SIFgpij;t+SlFusi;+SIFpaij = 
-[oa;] = —[V; - a]. (3.13) 
Points 4 


The above equation appears complex, however, it is merely an extension of the blade-only 
formulation (see equation 3.9). 

Hub and duct surfaces at the blade endpoints are streamlines of the flow, where the compo- 
nent of Viotai normal to the blade surface is zero. The hub and duct source panel influences on 
the blade vortex-lattice system induces zero CMV normal to the blade surface. This has been 
exhibited by a field-point velocity calculation routine at the control points on the blade surface 


[6], as well as by the velocity streamline traces developed in the next chapter. 
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Chapter 4 


Validation 


In this chapter, two methods of evaluating stators are compared. First, an infinite-pitch, con- 
stant 4 stator with a hub and duct in a uniform inflow is evaluated. The original PBD-14 
scheme will be compared to the source panel method for a variety of thickness distributions at 
different angles of attack to exhibit the differences in how each method represents of source- 
induced effects. Second, an infinite-pitch, constant $ stator in a uniform inflow is evaluated 
at the hub and duct to ensure their effects are accounted for properly. Finally, the number of 
blades will be increased to eleven and the effects analyzed at the hub. Since the most signifi- 
cant source-induced factors occur at this point, then cascade effect analysis will be performed 
here. The source panel method will exhibit its robustness when it uses varying blade thickness 


distributions for two different blade section geometries. 


4.1 The Infinite-Pitch, Constant 4 Stator 


This test case shows that the source panel method is more robust than the linearized slender- 
body theory of the original PBD-14 code for a larger range of thickness/chord ratios and number 
of blades. Results vary with vortex-lattice grid density and therefore all comparisons will be 
made with a 15 x 15 grid. Table 4.1 shows how the parameters for the stator will be varied in 


this test case. 
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Blade Section Shape | Angle of Attack (a) thickness/chord 


NACA66 modified 0° 26.67% 
Biconvex 12 20.00% 
13.33% 


2.67% 


Table 4.1: Infinite-pitch, constant t/D stator parameters 


Table 4.2: Infinite-pitch, constant t/D stator with a zero degree angle of attack. 


4.1.1 Initial Comparison for Large 5 Values at a Zero Angle of Attack 


Comparison of current PBD-14 and source panel method for the blade section shapes listed 
in Table 4.1 is exhibited in Figures 4-1 and 4-2. This set of figures shows how the linearized 
slender-body theory approximations for blade thickness distribution are inadequate when the 
blade section becomes thicker. Figures 4-3 and 4-4 show the results for the source panel method 
which has less than one percent error associated with the blade thickness distribution. Because 
of the straightforward shape (no camber, as seen in Table 4.2) and uniform inflow, this set of 
stators does not produce any circulation or forces which was verified by the PBDOUT.SGR file. 
The blade lattice structure is a 15 X 15 grid with hub and duct images (7) implemented. The 
dividing streamline is shown so that the modeled thickness distribution (velocity streamlines) 
can be compared with the actual thickness distribution on the pressure and suction side surfaces 


of the blade. The cut is made midspan on the blade. 


4.1.2 Comparisons at a One Degree Angle of Attack 


The next step was to place the blades at an angle of attack greater than zero. This would 
test equation 3.8 results that predicted the tangential velocity component on the mean camber 
surface (Ujn cosa) is the only affect the blade thickness distribution. Figures 4-5 through 4-6 
exhibit the continued trend for the original PBD-14 code to under predict the blade thickness 


distribution at the larger thickness/chord ratios. Also, the results of equation 3.8 are verified in 


43 


Figure 4-1: Linearized slender-body theory velocity streamline results for a NACA66 modified 
blade section with 4 = 0.020 and a = 0° in a uniform inflow. 


y 
0.20000000 
0.19999810 
0.19999239 
0.19998287 


0.00000000 
0.04999239 
0.09998477 
0.14997715 


0.00000000 
0.00087262 
0.00174524 
0.00261786 


Table 4.3: Infinite-pitch, constant t/D stator with an one degree angle of attack. 


Figures 4-7 and 4-8. Because of the angle of attack (Table 4.3 shows the stator’s B-spline net 
design for this a) and uniform inflow, this set of stators produces some circulation and forces 


which was verified by the PBDOUT.SGR file. 


4.1.3 Comparisons at a Five Degree Angle of Attack 


The next step was to place the blades at an angle of attack greater than one. This would test 
equation 3.8 results further. Figures 4-9 and 4-10 exhibit the continued trend for the original 
PBD-14 code to under predict the blade thickness distribution at the larger thickness /chord 
ratios. Also, the results of equation 3.8 are well supported in Figures 4-11 and 4-12, even though 


the streamlines do not exactly match the actual thickness distribution shape. However, this 
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Figure 4-2: Linearized slender-body theory velocity streamline results for a Biconvex blade 
section with — = 0.020 and a = 0° in a uniform inflow. 


y 
0.00000000 | 0.20000000 | 0.00000000 
0.04980973 | 0.19985742 | 0.00435779 


0.09961947 | 0.19971046 | 0.00871557 
0.14942920 | 0.19957226 | 0.01307336 


Table 4.4: Infinite-pitch, constant t/D stator with a five degree angle of attack. 


series of figures definitely shows the improvements that the source panel method makes to the 
overall blade thickness distribution representation in many different scenarios. Because of the 
angle of attack (Table 4.4 shows the stator’s B-spline net design for this a) and uniform inflow, 
this set of stators produces some circulation and forces which was verified by the PBDOUT.SGR. 
file. 


4.1.4 Varying the Thickness/Chord Ratio at a Zero Angle of Attack 


To show the robustness of the source panel method over the linearized slender-body theory, 
different thickness/chord ratios are compared in Figures 4-13 through 4-16, including the Figures 
4-2 and 4-4 (thickness/chord ratio of 26.67%). As as the thickness/chord ratio approaches zero, 
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-0.04 
-0.05 
-0.06 


-0.05 0 0.05 0.1 0.15 


Figure 4-3: Source panel method velocity streamline results for a NACA66 modified blade 
section with 4 = 0.020 and a = 0° in a uniform inflow. 


the source panel method and linearized slender-body theory begin to merge to the same blade 
thickness distribution, as expected. Note that as the blade thickness distribution approaches 
zero, it is better to increase the grid density to obtain a finer source and vortex strength 
distribution over the blade. In this manner, the blade representation is more accurate, but 


computational time increases. For some cases it may be necessary to increase the grid density 


to 25 x 25 or even 35 x 35. 


4.1.5 Hub and Duct Interfaces and the Cascade Effect 


Finally, the next example is just a check that the source panel method with hub and duct 
modeling yields the correct blade thickness distribution at these intersecting points on the 
model. Figures 4-17 and 4-18 exhibit the robustness of the source panel method to include 
all effects (blade-to-blade, hub-to-blade, and duct-to-blade) of the sourced-induced functions 
accurately. A special case of the hub interface with the blades is shown in Figure 4-19 where 
the number of blades increases to eleven. The source distribution increases to appropriately 


represent the cascade effect in the model. Note that the source panel method still exhibits 
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Figure 4-4: Source panel method velocity streamline results for a Biconvex blade section with 


4 = 0.020 and a = 0° in a uniform inflow. 


its robust nature to represent the blade thickness distribution at the hub and blade interface 


accurately. 
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Figure 4-5: Linearized slender-body theory velocity streamline results for a NACA66 modified 
blade section with 4 = 0.020 and a = 1° in a uniform inflow. 


Figure 4-6: Linearized slender-body theory velocity streamline results for a NACA66 modified 
blade section with 4 = 0.020 and a = 1° in a uniform inflow. 
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Figure 4-7: Source panel method velocity streamline results for a NACA66 modified blade 
section with 4 = 0.020 and a = 1° in a uniform inflow. 


Figure 4-8: Source panel method velocity streamline results for a Biconvex blade section with 


4 = 0.020 and a = 1° in a uniform inflow. 
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Figure 4-9: Linearized slender-body theory velocity streamline results for a NACA66 modified 
blade section with 4 = 0.020 and a = 5° in a uniform inflow. 
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Figure 4-10: Linearized slender-body theory velocity streamline results for a Biconvex blade 
section with 5 = 0.020 and a = 5° in a uniform inflow. 
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Figure 4-11: Source panel method velocity streamline results for a NACA66 modified blade 
section with $ = 0.020 and a = 5° in a uniform inflow. 


Figure 4-12: Source panel method velocity streamline results for a Biconvex blade section with 


4 = 0.020 and a = 5° in a uniform inflow. 
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Figure 4-13: Linearized slender-body theory velocity streamline results for a Biconvex blade 
section with 4 = 0.010 and a@ = 0° in a uniform inflow. 


Figure 4-14: Source panel method velocity streamline results for a Biconvex blade section with 


4 = 0.010 and a = 0° in a uniform inflow. 
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Figure 4-15: Linearized slender-body theory velocity streamline results for a Biconvex blade 
section with 4 = 0.002 and a = 0° in a uniform inflow. 
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Figure 4-16: Source panel method velocity streamline results for a Biconvex blade section with 


+ = 0.002 and o = 0° in a uniform inflow. 
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Figure 4-17: Source panel method velocity streamline results for a Biconvex blade section with 


+ = 0.015 and a = 0° in a uniform inflow at the hub/blade interface. 


Figure 4-18: Source panel method velocity streamline results for a Biconvex blade section with 


4, = 0.015 and a = 0° in a uniform inflow at the duct/blade interface. 
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Figure 4-19: Source panel method velocity streamline results for a Biconvex blade section with 
4 = 0.015 and a = 0° in a uniform inflow at the hub/blade interface. This streamtrace is for 


the key blade of an 11-blade stator system. 
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Chapter 5 


Conclusion 


5.1 Source Panel Method Results 


In this thesis, several enhancements to the blade thickness representation for lifting-surface 
design and analysis for marine propulsors were discussed. Optimization over the whole propulsor 


allowed for: 


e more accurate representation of the blade thickness distribution. This was accomplished 
by incorporating the effects of blade-to-blade interaction (blockage effect), as well as the 
hub-to-blade and duct-to-blade influence effects. 


e still able to use the original PBD-14 code structure to obtain solutions since the blade 
thickness distribution is represented by constant-strength source panels whose effect is in- 
corporated into an overall source-induced velocity term on the right-hand side of equation 


2.6. 


e more accurate results for larger ranges of thickness/chord ratios and numbers of blades 


(cascade effect). 


The weakness of the source panel method is that the computational time increases when a 
solution of linear equations is sought to determine the source strength distribution. Also, as a 
finer grid density is used to obtain more accurate results, the computational time increases. 

However, this panel method enabled the kinematic boundary condition to be satisfied on 


the actual pressure and suction side surfaces of the blade. This ensures that Vi,:a; on the blade 
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surface is tangential. The hub and duct source influences, as well as other blades influences, 


were shown in the validation results of Chapter 4 for the following conditions: 


| e varying thickness/chord ratios, 
e vary angles of attack and 


® two distinct blade section geometries. 


The velocity streamline traces showed that the source panel method is more accurate and 
robust over a the specified range of thickness/chord ratios, angles of attack and blade section 


geometries. 


5.2 Recommendations for Future Enhancements 


This thesis implemented a panel method into the PBD-14 code along the mean camber surface 
of the vortex-lattice structure. Some recommendations for further improvements and future 


work in this area include: 


e integrate the source panel method with the vortex lattice hub and duct formulation devel- 
oped in PBD-14.4 coding. This would allow for a more accurate representation of the hub 


and duct, as well as represent their associated affects on the blade thickness distribution. 


e develop a more accurate field point velocity grid calculation for the entire propulsor vice 


just a chordwise strip to capture all three-dimensional effects visually. 


e carry the source-induced velocity vector from the panel method solution to calculate a 
more accurate force and pressure distribution over the actual pressure and suction side 


surfaces of the blade. 


ov 
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Appendix A 


Constant-Strength Quadrilateral 


Source Distribution Formulation 


fi 


Rectilinear Panetf if (a ~ TO 
t ip ie ele ‘ \ 


Lf 
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P 3 
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o(%y) = constant 


(X15 Vir 8) 


Figure A-1: Rectilinear to Quadrilateral constant-strength source panel transformation. 


Consider a surface element (z = constant) with a constant-strength source distribution, 04, per 


area, bounded by four straight lines as described in Figure A-1. The potential at an arbitrary 


point P(z,y,z), due to this element is 


oA dS 
% __oA dS Al 
on) =-2 SS Sa 


and the velocity components can be obtained by differentiating the velocity potential: 


Ob OD AG 
(u,v, w) se (an Oy’ De? (A.2) 
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A.1 Near Field Formulation 


Execution of the integration within the area bounded by the four straight lines requires a lengthy 
process. The details of integration are shown in reference [12] and the results obtained for the 


potential of a planar element becomes 


(z-21)(y—y1) wv yi)(z2 #1) yp Tautratdie 


ritre—di2 
aos 22)(y3—- nt y2)(23— 2) In ratr3+da3 
rotr3—da3 
ae 23)(ya mi we y3)(@a—@3) y Tatratdaa 
dg rg+ra—d3a 
(z— za)yi-— valu ya) (21 — 4) | rat+ritda1 
CA at DD ytr—day 
O(x,y,2) aaa ie as (A.3) 
4a tan —1 f mizei—hy —tan7! my2eg—hg 
271 27Q 


+tan7! (macz=le ) acpane l (mats | 


+ |2| : 
+tan7! (masgaphs ia | —tan7! ( msscazhe ea ) 
+tan7! (magante ) — tan! ( maeimt ik ) 
where 
diz = 4/ (x2 — 21)? + (y2-— 1)? (A.4) 
doz = 1/ (x3 — x2)? + (y3 — yo)? (A.5) 
dg4 = 4/ (aa — %3)? + (ya — ys)? (A.6) 
day = 4/ (21 — v4)? + (yi — ya)? (A.7) 
and 
my = 2 (A.8) 
2 — XH 
m3 = 32 (A.9) 
3 — XL 
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Mm34 = vas (A.10) 
L4 — 3 


May = ous ue eR (A.11) 
LB — Lh 
and 
Th = 4/(%— te)? + (y — ye)? + 2? (A.12) 
ey = (a — xy)? +2? (A.13) 
hy = (@ — Le)(y — Yr) (A.14) 


where k = 1, 2,3, 4. 
The velocity components, based on the results from Hess and Smith [12] are 


Y2=W Jy 7 +re—di2 + Y8—¥2 Jy rg+r3—da3 
tA dio rytro+die d23 rgtr3+d23 (A.15) 


An 4 WA=U8 Jy Tatramdag 4 Yi—¥4 Jy Tat ri —day 
d34 rg+ratd34 dai r4a+rit+dai 


E1—22 Jy 71 -+re—di2 sie 22-23 Jy ratr3—do3 
v= oA dia ritret+di2 d23 To+r3+da3 (A 16) 


Ar 4 23—%4 In rgtra—dza a Za-Z) Jy Ta +rj3—da4} 
dg rg+ra+dg4 day ratritdai 


tan} (muerte 2€ = ) —tan7! (magenta ea—h ) +tan—! (maseamla ) —tan7! (maaga=ts esha ) 


= OA 2ry 2r2 Zr3 
4n ae tan—! (masee='s ) = tan! (masge=te ) a tan! ( sents ) _ tan (arent ) 
(A.17) 


The w and v components of the velocity are defined everywhere but at the edges of the 
quadrilateral, where they become infinite. If the influence of the element on itself is sought, 
then these velocity components approach zero near the centroid. The jump in the normal 


velocity component as z —> 0 inside the quadrilateral is 
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w(z = 0+) = +4. (A.18) 


When the point on interest P (z = constant) lies outside of the quadrilateral, then 


we 02) = 0; (A.19) 


A subroutine was written to incorporate code to handle nonplanar panels. Basically any rec- 
tilinear panel is input into this subroutine which then adjusts the panel dimensions to correctly 
represent it as a quadrilateral panel with a x —y planar surface on a constant z— surface. This 
subroutine also calculates the center of the rectilinear panel which maintains the orientation 


between the local (panel) and global (blade geometry) coordinate systems. 


A.2 Far Field Formulation 


For improved computational efficiency, when the point of interest P is far from the center 
of the element (29, yo,z = constant) then the influence of the quadrilateral element with an 
area of A can be approximated by a point source. When the term “far” is invoked, that is a 
distance more than five average panel characteristic lengths (longest distance between any two 
endpoints), then a simplified approximation is used. The point source influence for the velocity 
potential becomes 

oAA 


A (CET OT a — 


The velocity components of this source element are 


oAA(x — Xo) 


i Flea eG ee (A.21) 
Sp AA) 2 cs 
An [(x — 20)? + (y — yo)? + 23/2 (A.22) 
= rai (A.23) 


4n [(x — 20)? + (y — yo)? + 22)°” 
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An algorithm for calculating the influence of a quadrilateral constant strength source element 


(Program 11) was used [11]. 
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Appendix B 


Infinite-Pitch Stator PBD-14 Input 
File | 


PBD-14.3 admin file for Infinite-pitch stator DSH 6/5/00 


stator. bsn ‘File name for blade b-spline net 
samp.vel :File name for wake field 
1 15 15 :nblade, nkey, mkey 
1 1 rispn (O=uniform,1=cos),iffixlat 
14 1,2, .., 40 :mctrp,mc(n) 
7 0.0 7 0.0 tihub, hgap, iduct, dgap 
0 0 7IHUB_SUBLAY, IDUC_SUBLAY 
10 8 ‘HDWAK, NTWAKE 
0.84 300 :Cq, MXITER 
0.3 0.3 :OVHANG(1), OVHANG(3) 
10—-10 0 8 ‘nx, ngcoeff, mltype, mthick 
4 -IMODE 
0 :NXIMAX 
5 001 .02 5. 1 :niter, tweak, bulge, radwegt, nufix 
1.08 snplot, hubshk 
4 2 :NOPT, NBLK 
999. 10 15  .05 :ADVCO, XULT, XFINAL, DTPROP 


o[ o| 0, 0] 0, 0] Oo, oy] 0, oT 
27a] 4, 5, 6| 7, 8] 9| 05] 10] afk 
1/D 
of of of of of of of of of o[ca 


Table B.1: Sample input file for PBD-14. 
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